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Abstract 



The nuclear sigma term is calculated including the nuclear matrix element of 
the derivative of the NN interaction with respect to the quark mass, m q 9 q^ n ■ 
The NN potential is evaluated in the skyrmion-skyrmion picture within the 
quantized product ansatz. The contribution of the NN potential to the nuclear 
sigma term provides repulsion to the pion-nucleus interaction. The strength 
of the s-wave pion-nucleus optical potential is estimated including such con- 
tribution. The results are consistent with the analysis of the experimental 
data. 
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I. INTRODUCTION 



Nowadays Quantum Chromodynamics is accepted as the theory for strongly interacting 
particles. It was probed with success at high energies where perturbative expansion can 
be used. Unfortunately in the low energy regime the perturbative expansion is not valid 
and nonperturbative schemes such as lattice QCD have been implemented, but still with 
enormous computational efforts. The energy region in which quarks and gluons build the 
effective hadronic degrees of freedom is of interest for nuclear physics applications. Effective 
theories have been constructed considering only hadronic degrees of freedom and designed to 
work at low energies. In such models the main guides are the symmetries of the fundamental 
theory. A bridge can be made between the two descriptions when considering the QCD 
vacuum structure. The presence of hadronic matter affects the vacuum structure, which is 
reflected by the modification of the quark and gluon condensates PJ. The hadronic matter 
dig a hole in the local value of the condensate as expressed by the expectation value of the 
local composite operator qq, where q is the quark operator. 

The difference between the volume integral of the quark condensate in the presence of 
hadronic matter and its vacuum value is proportional to the nuclear sigma term. It can 
be obtained by the application of the Hellmann-Feynman theorem, which allows to express 
the nuclear sigma term as mq^f- with Ma being the nuclear mass We are limiting 
our discussion to the SU(2) flavor sector. The nucleon-nucleon interaction is responsible for 
the nuclear binding, so it is natural to question how the NN potential will be affected by 
variations of the quark masses, if a reliable calculation of the nuclear sigma term is desired 



However, to predict how the potential depends on the current quark masses of the fun- 
damental theory, the connection of the effective hadronic model with QCD is necessary. 
The quark mass terms break explicitly the chiral symmetry in the QCD Lagrangian. In 
the SU(2) flavor sector of the theory, the up-down quark masses originate the pion mass, 
since the pion is considered the Goldstone boson arising from the spontaneous breaking of 
SU(2) L xSU(2) R chiral symmetry §. 

The Skyrme model allows to study the effect of the explicitly chiral symmetry breaking 
in a chiral invariant hadronic model. The quark mass comes through the pion mass term in 
the Skyrme Lagrangian. Such connection is made using the Gell-Mann, Oakes and Renner 



where m q = (m u +md)/2, m u and are the up and down current quark masses, respectively. 
The mass m q is quite small in the hadronic scale and the local value of the vacuum condensate 



The Skyrme Model describes qualitatively the nucleon || and the nucleon-nucleon 
interaction fTjHH . The particular feature of the model is that the pion mass term in the 
Lagrangian contributes not only for the nucleon properties but also to the NN potential. 
We are supposing that the dependences of the pion weak decay constant (/„■) and the Skyrme 
parameter (e) on the up-down quark masses are small and they are not considered here. 

The derivative of any hadronic quantity, O, with the up-down average quark mass in the 
Skyrme model is translated into a derivative in m^. using the GOR relation, Eq. ([!]), such 
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that 



80 2 dO 

m„— — = m — — -. (2) 

om q om^ 

In this work, we obtain the derivative of the NN interaction with respect to the cur- 
rent quark mass by calculating the potential with the product ansatz approximation in the 
Skyrme model. We discuss how the derivative enters in the evaluation of the nu- 

clear sigma term in the context of nonrelativistic quantum mechanics and its contribution 
to the quark condensate. The expectation value of in the nuclear state contributes 

to the repulsive strength of the low-energy isoscalar pion-nucleus s-wave optical potential. 
These calculations can be tested through the comparison with the results of the pion-nucleus 
scattering and pionic atom data analysis fTTl|-|T2|| . For completeness, we also calculate the 
nucleon sigma term, the derivatives of the pion-nucleon coupling const ant (g^NN) and axial 
charge (g A ). 

In Section II, we discuss the nuclear sigma term in the context of nonrelativistic quantum 
mechanics, defining the part of the sigma term associated with the NN interaction. In Section 
III, we discuss the quark condensate at finite nuclear density with the inclusion of the new 
term. In Section IV, the Skyrme model is briefly discussed. In Section V, the results for 
the derivative of the NN potential, with respect to the quark mass, are presented together 
with the derivatives of the one-body observables. In Section VI, part of the main results 
of this work are presented, those correspond to the calculations of the nuclear sigma term 
for several nuclei. In Section VII, the contribution of the derivative of the NN potential to 
the pion s-wave optical potential is calculated and compared to the results of the analysis 
of experimental data. In Section VIII, we summarize our conclusions. 



II. THE NUCLEAR SIGMA TERM. 



The sigma term is a measure of the isoscalar content of the nuclear state {^ A ), with A 
nucleons, and thus it gives the strength of chiral symmetry breaking due to the u and d 
quark masses. It can be directly defined as the expectation value of the commutators of the 
axial charge (Q 5 ) with the Hamiltonian [13,14], given by 



Q\ Q 5 ,H 



In QCD framework it can be expressed as [|14 



2m g / d 3 x( {^ A \q(x)q(x)\ ^ A ) ~ (0\q{x)q(x)\0)) = m, 



dM A 



(3) 



(4) 



it is proportional to the nonstrange quark condensate. The last equality in Eq. ([|), follows 
from the Hellmann-Feynman theorem [|l[] at the level of the fundamental QCD theory, where 
the quark masses contribute to the Hamiltonian as SHqcd = m u uu + m^dd. 
The nucleon sigma term is 



Eat = m q 



dm„ 



N 



(5) 
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where we used (0) and Mjy is the nucleon mass. The experimental value extracted from 



pion-nucleon scattering is Ejy ~ 45 MeV |15 



The nuclear sigma term has a contribution from the derivative of the nuclear binding 
energy {E B ) with respect to the quark mass. Using that M A = ZM P + {A — Z)M n + E B , 
where Z is the proton number, M p and M n are the proton and neutron masses, respectively, 
we have: 

S A = AY, N + m v^f- ' ( 6 ) 

where the neutron and proton sigma terms were considered equal. 

The effect of the nuclear binding in the sigma term can be calculated, for example, from 
the Schrodinger equation in the center of mass system. If the quark masses contribute 
effectively to V^n, they also influence the binding energies in the Schrodinger equation, 
written as 



H\* A ) 



i=1A ZM N I ^. 



\* A ) = E B \* A ) . (7) 



Using the Hellman-Feynman theorem and considering that Mn and V%j depend on m q , 
the nuclear sigma term is given by 

S A = AL N + Ef + TX A , (8) 

where 

is the contribution of the NN potential to the nuclear sigma term, and 

Si = ~^Ek , (10) 

M N 

is the term coming form the kinetic energy Ek, which is the expectation value of the sum 
of the nucleon kinetic energies in the center of mass frame of the nucleus. The AH n term 
in Eq.(|8|) has been discussed in the context of nuclear matter in Ref. The last term in 

Eq. (H) can be rewritten using Eq.(|2]), as 

1^ dVij l v 2 dVij 

- / m„- — - = - > mi— — - . (11) 
2^ q 8m q 2^ dml V ; 

The above equality permits to calculate the derivative m q ^f - in the Skyrme model. 

III. THE QUARK CONDENSATE 

The nuclear sigma term is responsible for the change in the quark condensate at finite 
density [3]||. The sigma term for nuclear matter at a finite density (p) in a given volume 
{Vol) is given by 
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E p = 2m q Vol( (p\q(x)q(x)\p) - (0\q(x)q(x)\0)) = m q ^^- , (12) 

where £ is the nuclear matter energy density given by £ = p + Ev+Ek p, in the limit of 
A — ► oo and Ey is the potential energy. Since p = A/Vol, we get from (O) 



Using (D, the definitions E^©, ^a© and the GOR relation, 2m q {qq) vac = -m 2 Jl, we 
have 

~ 1 ' p + ... . (14) 



If we neglect (t,\ + E^- ) the expression above reduces to the leading order term found in 
Ref. ||. Recent works by Li and Ko ||16|| , Delfino et al 0] and Brockman and Weise Jl8 



calculated the quark condensate in relativistic nuclear matter theories following the work 
of Cohen and collaborators ||. They included in the evaluation of the nuclear sigma term 
the derivatives of the meson masses with respect to the quark mass, and in their context 
the quark condensate didn't deviate strongly from the leading order calculation at normal 
density. Brown and Rho [[19] have studied the quark condensate in nuclear medium without 
the contribution of the NN potential to the nuclear sigma term. The values of the condensate 
at normal density, in all cases, are around the leading order result, (qq) PN / (<?<?) mc ~ 0.70 . 

In the following section we will calculate E^ + E^ and E^ using the Skyrme Model and 
the GOR relation. Subsequently we will calculate the contribution of to the pion-nucleus 
optical potential. In the summary, we will present an estimate to the quark condensate using 
(p~4]), with E^ obtained from the Skyrme model. 



IV. THE SKYRME MODEL. 

The Skyrme Lagrangian M can be expressed as 

+ l -flml{TrU-2) . (15) 

where U = e tT ' 7T ; f are the Pauli isospin matrices. Assuming the hedgehog ansatz, n (r) = 
fF(r), and minimizing the energy we arrive at a differential equation for the profile function 

(F) ' ' ' ' 

(^. + 2s 2 ) F" + -uF> + F' 2 2sc - ~sc - ^ - /?VsinF = , (16) 
V 4 J 2 2 u 2 

where u = 2ef n r, [3 = m n /(2ef n ), s = sinF and c = cosF . The solution with baryon 
number B=l, can be achieved imposing F(0) = n and F(oo) = 0. 

We make use of the quantized version ]7|,|8| of the product ansatz to get a qualitative 
idea of the consequences of finite quark masses in the NN potential. The product ansatz 



J 7T / 



Tr (d„Ud»U ] 



-Tr 



rfdJJ, U^d u U 
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was originally introduced by Skyrme || as an attempted solution for the skyrmion-skyrmion 
interaction. It has the virtue of giving baryon number B=2 for any separation distance and 
gives the one pion exchange potential for large distances. It provides an NN potential which 
can be decomposed in the form 

Vnn = V c + fi.f 2 [Vss (01-02) + V T S 12 ] , (17) 

where Vc is the isoscalar central potential, Vss is the spin-spin potential and Vp is the tensor 
potential; S12 = 3 (a\. f) (oV^) — aV^V ; a are the Pauli spin matrices. 

Exact numerical solutions for the skyrmion-skyrmion interaction were obtained in the 
work of Walhout and Wambach |§. In that work ||, they compared the exact numerical 
calculations to the product ansatz result and to the Paris potential. Although the gross 
feature of the central component is given by the product ansatz, it misses the attraction 
which is present in the exact calculation and Paris potential. The product ansatz is roughly 
reliable for r > 1.5 fm in the tensor and spin-spin channels. 

Although the product ansatz in the NN central potential is only approximated when 
compared to exact numerical calculations Q, we believe that is possible to use it in the 
calculation of the derivative. According to Jackson et al 0, the difficulties in reproducing 
the NN interaction from the product ansatz come from the fact that the potential is a 
difference of two large quantities of magnitudes two orders higher than the potential itself, 
and so little deformation in the skyrmions can originate drastic effects. This shall not occur 
in the derivative of the potential with respect to the quark mass, since the final result is of 
the same order of magnitude of the individual nucleon contributions. So, deformations that 
could occur in the skyrmions as they are close, would be in this case only corrections that 
could not change the magnitude of the result. 



V. DERIVATIVES OF THE NUCLEON OBSERVABLES AND NN POTENTIAL 

We begin by showing the results of the derivatives of the one-body observables with 
respect to the quark mass. We can calculate several nucleon observables 0, such as the 
nucleon mass, the pion-nucleon coupling {q^nn) and the axial coupling (<m), from the Skyrme 
Model in the semi-classical quantization approach. Thus, by changing the pion mass we can 
easily calculate m^^f, m% d g^2 N and J^|- • These calculations provide the determination 

of the sigma term, m q 99 q^ n and m q ^^, respectively. The results of the derivatives of the 
observables are presented in Table I, and for completeness, the value of the observables 
themselves are also shown. 

The calculations were performed with two sets of parameters, /„- = 54MeV, e = 4.84 || 



and /„- = 93MeV, e = 4.0 [20] with m n = 138MeV. In Table I, the values of the nucleon 
sigma term, were included. It has been calculated previously in Ref. ||. They obtained 38 
MeV and 49 MeV, for the set /„- = 54 MeV and e = 4.84, using different methods without 
recurring to the Hellmann-Feynman theorem. These values are comparable to our result of 
52.2 MeV. The experimental value of about 45 MeV ||15|| , is somewhat consistent with our 
value of 59.6 MeV for /„. = 93 MeV and e = 4.0 (see Table I). We will restrict our evaluation 
of the nuclear sigma term to this set of parameters, since /„- is taken from the experimental 
value. 
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The coupling constant g W NN has very different derivatives for the two sets of parameters, 
while the derivatives of qa are similar, as shown in Table I. However, in the case of the 
9-kNNi the derivatives m q d§ Q^ N , the nucleon sigma term and are consistent with the 



Goldberger-Treiman relation [|4 

The derivative is evaluated numerically in the Skyrme model using the quantized 

product ansatz with the parameter set f n = 93 MeV and e = 4.0. We calculated the NN 
interaction in skyrmion-skyrmion picture and recalculated it for a small change in the pion 
mass. Thus, we obtained numerically the quantities m 9 |^-, m q §^- and m q ^^- • They are 
shown in Figs. 1, 2 and 3, respectively. The spin-spin and tensor potential derivatives are 
compared to those extracted from the one pion exchange potential (OPEP), presented in 
the appendix, where qa and rn q ^- were taken from Table I. 

The results for function of the relative distance of the skyrmion centers, 

are shown in Fig. 1. The calculations show minimum around 1 fm, with values about -20 
MeV, while at short distance they present repulsion. The difference y( ni7 '- Vi ^ _ y^™ 1 *- ) j g 



comparable to m g |^, as shown in the figure. 



r)V dV OPEP 

In Fig. 2, the results for rriq-Q^- are shown. The derivative m q — S- — is compared with the 
corresponding one for the skyrmion-skyrmion interaction. They agree for separations of the 
nucleons above 3 fm. This is expected since the OPEP tail is present in the skyrmion- 
skyrmion interaction. The derivative found with OPEP is negative at short-distances, 
while the product ansatz results are positive. The shapes of the curves for the difference 
_ yim^-o) m ^^_ are qualitatively similar, although the tail presents a no- 
ticeable difference. We have checked that at distances above 3 fm, the difference in the 
potentials with finite and zero pion masses is given by the results from OPEP, since the tail 
is dominated by being zero or finite. 

The derivative of the spin-spin component of the skyrmion-skyrmion interaction, m q ^^-, 
is shown in Fig. 3. The derivative of the OPEP spin-spin component for distances above 3 
fm agrees with the present calculation. As expected for distances below 2 fm, the derivatives 
of the skyrmion-skyrmion potential and OPEP become quite different. The shape of the 
difference Vg™* - — vj™'" -0 '' is similar to that of ™ q ^f- ■ However at distances above 3 
fm, the difference in the potentials with finite and zero pion masses is essentially dominated 
by the results from OPEP, which explains the origin in the differences found. 



VI. RESULTS FOR THE NUCLEAR SIGMA TERM 

One way to extract the nuclear sigma term is from the low energy pion-nucleus scattering 



analysis. In fact, in a recent work of Brown and Rho |19j it is suggested the possibility of 
obtaining information about the quark condensate from pionic atoms. 

The NN potential extracted from the product ansatz is not realistic in the sense that it 
does not reproduce the central attraction in the mid-range. This cannot be ignored since it 
means that there is no wave function for the bound state in this potential. Then we decided 
to give up consistency and from the derivative of the calculated potential, we can estimate 
its contribution to the nuclear sigma term using charge densities from the electron elastic 



scattering 21 
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In the particular case of the deuteron, we calculated 



dV 

S^H = (^ 2 H| m <?^~ l^ 2 H) ( 18 ) 



using the normalized wave function f22|] . 



= = ^y 1 ' 2 + ^^ 12 (r)} XlM. (19) 

where w and w are the s-wave ( 3 Si) and d-wave ( 3 Dx) from the deuteron wave func- 



tion obtained with the Reid potential [23]. The normalized wave function is such that 
Jo°° dr (u 2 + w 2 ) = 1. The spin component Xim = \1M) is such that Ximx\m = 1 an d 
r = |ri — r 2 | is the relative distance between the nucleons. 

The Skyrme model offers naturally forces of two, three and more bodies. Nevertheless, 
it is still a good approximation to neglect them all in the nuclear sigma term except the 
contributions from the two-body potential. Assuming a modified Hartree approximation 
such that it has the correct number of pairs, we have 

s y = ma^ j dhip(ri) j^ ^-i) _ (20) 

where / <i 3 rp(r) = 1 and A is the mass number. 

For isoscalar nuclei heavier than the deuteron we neglected the spin-spin and tensor 
components of the derivative of the potential, considering that they are more than one order 
of magnitude smaller than the central part and the isospin is averaged in the matrix element 
(which yields zero in the Hartree approximation). Then only the derivative of the central 
potential enters in the estimate of H\. 

We present in Table II the results of E^ using ( p0|) and the densities obtained from 
elastic electron scattering f2~I| . In this first approach, we assumed the same proton and 
neutron density distributions. The deuteron results were included in Table II. In the case 
of the deuteron the contribution from the NN interaction is — 3.60MeV . We expected that 
we could roughly count the number of NN pairs for light nuclei in E^. For instance, in 4 He 
we can count 6 pairs, then we should have 6x (— 3.60)MeV=— 21MeV that is a lower value 
(in modulus) than the calculated — 34.5MeV. A possible explanation is that the deuteron is 
"big" and its NN interaction is on an average weaker than in 4 He. Triton has 3 pairs and 
we estimate the contribution of the potential to the nuclear sigma term as being between 
3xE^ H = 10.8MeV and 3xSr Hc /6 = 17.2MeV. Averaging we find a value of 14MeV for 
triton. The same should work for 3 He. For larger values of A, remains always smaller 
than the number of pairsxE^ H or the number of pairs xEl^/6, because the NN interaction 
has a range around 1 fm. Also we observe that E^/ 'A is always growing with A except for 
the 197 Au and 208 Pb. This can be explained by considering that lead have lower density 
distribution at mid-range than gold, as we checked. 

In Table II, we also included the kinetic energy contribution to the sigma term. The 
kinetic energies from the deuteron, 3 He and 4 He were extracted from the work of Schiavilla 



et al [24] . In the heavier nuclei we took the average kinetic energy of nonrelativistic nucleons 

4i 



in nuclear matter, given by Ek/A = j^(2)ej. At p = 0.17fm 3 , kf = 1.36fm 1 ,Sf =38MeV 
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and we arrive at Ek/A = 22.8MeV. The values of Tj A + = + are also presented and 
they will be used in the evaluation of pion-nucleus optical potential. 

We noticed that the calculation of the sigma term from the NN potential using the 
Gell-Mann, Oakes and Renner relation(GOR) and taking the derivative with respect to 
the pion mass, produced very different results from those calculated considering only the 
symmetry breaking term (according to Adkins and Nappi ||). This occurs because there 
is an enormous contribution from the fourth order term. We also observed that it is in the 
mid-range (lfm-2fm) that the potential NN gives the biggest contribution to the sigma term. 



VII. THE NUCLEAR SIGMA TERM IN THE PION-NUCLEUS S-WAVE 

POTENTIAL 

Considering that in the last section we got an estimate of the contribution of the NN 
potential to the nuclear sigma term, we can also use it in the calculation of the s-wave 
pion-nucleus potential. The interaction Lagrangian between the pion and the nucleus which 
is originated by the nuclear sigma term(S^) is given by 



s 

int. 



m, 



^1 (ip A \ : q(x)q(x) : |^a)tt 2 . 



(21) 



We can approximate 



m g 0>a\ ■ q{x)q(x) : \ip A ) 



A 



'P(r) 



(22) 



with / d 3 rp(r) —A. The nuclear sigma term is separated in its contributions of individual 
nucleons, potential and kinetic energy; = AT, N + + (see Eq. H). 

The component of the interaction Lagrangian which originates the s-wave potential 
(U l=0 (r)) can be written as 



C 



1=0 

int 



n 



1=0 1 



-IT 



(23) 



where the s-wave self energy is given by 



n 



1=0, 



-Mi + m7T /M N M2 e + 



1 + 



4nf 2 A y M 



N 



'Mr)- 



(24) 



We have separated and b^ e contributions in the self energy Il i=0 (r), Eq. (pl[) . 

This last one should contain all the mechanisms of the pion-nucleus interaction not included 
in T, A r+K . The normalization was chosen such that bl™ e has the same normalization as the 
ttN scattering length. 

Introducing the self energy, Il z=0 ( 
uj 2 = k 2 + m 2 and dividing by 2m v , we get 

v 2 2uVn md V2„, n l=0 < 



r), in Klein-Gordon equation for the pion p2|,^5j, using 



2m T 



+ 



2m„ 



V 2 

v Coul 

2m^ 



+ 



2m„ 



<Pa(r) 



2m. 



(25) 
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In this way we have obtained the pion Schrodinger equation, where we identify the potential 



U l =°(r) = ^n'=°(r). (26) 
The part of the s-wave optical potential from the term is given by 

= Sa^ (27) 

where Eq.s (|4]) and ( p6|) were used. Then, if T^ +K is positive, it will produce an attractive 
contribution to the optical potential and repulsive otherwise. 

The scattering amplitude in the long-wave length limit and in the Born approximation 
is given by 

^ f j3„„ik.rrrZ=0/ 



/(k) = / d 3 re^ r U l=0 {r). (28) 

The nucleus recoil is considered in Eq.([28l) by the reduced mass factor /x n = m 7r M J 4/(m vr + 
Ma), that corresponds to the addition of the nonrelativistic kinetic energy of the target in the 
pion Schrodinger equation (p5|). The pion-nucleus scattering length in Born approximation 
is given by a = /(0). 

Although we have discussed the Born approximation in general, it only will be used 
for the deuteron, 3 He and 4 He, when discussing experimental data. For heavier nuclei 
we will compare the calculations to the Seki-Masutani parameter 6q extracted from 

experimental data analysis and defined as the coefficient of the p(r) term in the optical 
potential, 

2m^f/ /=0 (r) = -4?r(l + m n / 'M N )b s M p(r). (29) 
Comparing the equations (p^) , (p6|) and (|29|) , we get 

bl M = bC + bV K , (30) 



where, 



T v+K / m 

r = ^i+^ (3D 



IN 



which is the contribution of the nuclear sigma term from the NN potential plus the kinetic 
energy to 6q M . 

The isoscalar scattering length fe^; 6 was obtained using the experimental value of the ird 
scattering length (Rea^ = -0.0264(ll)m; 1 ) [|6|, 



<: = ( i 



M 



N 



v M~J~2 ~8^jf 



(32) 



with m n = 138MeV, M N = 938MeV, M d is the deuteron mass, f n = 93MeV. We used 
^Y^ K = — 4.46MeV, extracted from Table II, for f w = 93 MeV and e =4.0 . These parame- 



ters adjust the irN p-wave scattering f20fl . Substituting these values in ( p2|) we get 
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b f C = -0.0099 (m; 1 ). (33) 
This value is consistent with —0.0077(11) m" 1 got from Sigg et al |E7J| in the isospin 



symmetry hypothesis and also with the value —0.0083(38) m~ from Koch [28 

It has been noted that the scattering length for bound nuclei is shifted from the scattering 
length predicted by the theory for free nucleonsQ. This suggests that the cause for this shift is 
precisely the nuclear sigma term and kinetic energy contribution. The value of the isoscalar 
tcN scattering length parameter b is controversial and is not our purpose to calculate it here. 
Then, we obtained the effective isoscalar 6q^ £ originated from experimental pion-deuteron 
scattering length that contains the single and double scattering processes not included in 

spV+K 

We constructed Table III using the results of the potential and kinetic energy contribu- 
tions to the sigma term in Eq. (|30|) . The experimental results were extracted from pionic 
atoms and pion-nucleus scattering analysis []10"| , |IIf . 

In the particular case of 3 He, b^f 1 was evaluated subtracting the isovector contribution, 
i.e., considering only single scattering, 



u cB tv' A A 



(34) 



and using the experimental values Rea^- 3 He =0.056(6) |22[] and b\=— 0.0962(7) extracted 
from [27| we got fro"? ( 3 He)~ —0.015. The experimental result for 12 C was obtained by the 



Jeff 

16/ 



Seki-Masutani analysis and also for O extrapolated from tt scattering data at 29MeV and 
50 MeV QT0j| . The experimental result for 40 Ca was found in the Masutani-Seki analysis fll 



extrapolated from n scattering data at 25MeV to zero energy. In the cited work [|HJ, it is 
also found that the pion interacts with the nucleus at an effective density p e = px/2. The 



empirical value of — 0.0304m~ 1 was obtained by Salcedo et al |T2| averaging values obtained 
from the analysis of pionic atom and pion-nucleus scattering data. 

Our evaluation of the nuclear sigma term from the NN potential generates a repulsion 
in the pion- 12 C optical potential in the central region which is approximately U^ +K (r) = 

T, v+K 

~ 2m / 2 A ^( r ) ~ ^ MeV that is roughly of the same order necessary to fit the pionic atoms ( 
^15to 30 MeV) and is "weakly dependent of A" @ . 



VIII. CONCLUSIONS 

We have calculated the derivative of the skyrmion-skyrmion interaction with respect to 
the quark mass, m q 9 ^ N , in the quantized product ansatz. We have related it to m ^^f 



^From Kluge p0[] : "The conclusion is that this repulsion cannot be obtained by an iteration of 
irN interaction, but rather represents an independent feature of the interaction of pions with bulk 
nuclear matter. An attempt to explain this as a binding effect of the nucleons has been undertaken 
in the framework of relativistic mean field theory of the nucleus (Birbair et al 1983 J.Phys.G: 
Nucl.Phys 9 1473; Goudsmit et al 1989 Preprint ETHZ-IMP RP-90-03 ETH Zurich)" 
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using the Gell-Mann, Oakes and Renner relation, which made the calculation possible. We 
compared it to the OPEP calculation, and they agreed for distances above 3 fm. The matrix 
element of above derivative of V^n m the nuclear state is part of the nuclear sigma term. The 
magnitudes of the quark mass contribution to the nuclear potential were found appreciable. 
We calculated the contribution to the nuclear sigma term arising from the nuclear interaction, 
assuming that the nuclear density in finite nuclei is equal to the experimental charge density 
conveniently normalized. 

The contribution from the NN potential to the sigma term was incorporated in the effec- 
tive isoscalar scattering length of the pion-nucleus optical potential. There was a remarkable 
agreement with the empirical values extracted from the analysis of pionic atoms and scatter- 
ing experimental data. We showed that it is possible to explain the repulsion in the isoscalar 
channel of the pion-nucleus optical potential, if we take into account in the nuclear sigma 
term the contribution of the NN potential. 

In nuclear matter, we should include the contribution of the NN potential in the 
nuclear sigma term in the form = AH^ + + S^. This correction alters the value of 
the condensate (qq) p / {qq) vac = 1 — P^a/ Am^fl, and as is negative there is a tendency 
to a drastic reduction of XU- Then the condensate at normal nuclear density is about 
(<?<?) PJV / (Q^vac ~ 0.95, and by this way it does not decrease as much as it was calculated 
before. More definitive conclusions about the quark mass dependence of the NN potential, 
should be pursued in exact evaluations of the skyrmion-skyrmion interaction, that is beyond 
the present work. 
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Appendix 

The nonrelativistic reduction of the one pion exchange interaction is given by p9| 



yOPEP = 9a (rnA 2 ^ 



- (^.ay + ( — 2 + — 

3 \(m n r) m^r 




4tt \2f n/ 

V s °s PEP (^.a 2 ) + V T OPEP S 12 , (a.l) 



\Tl-T2 



The derivatives of the potentials with respect to the square of the pion mass can be 
simply stated as 



2 

dml 12n\2fJ \~^dml"^ ' * A V 2 



m n Q o ~ TcT~ I WF~ 1 S 2 9ATT^ m n + 9a U--S- ? ( a - 2 ) 



and 



2 

"* dml 4tt \2fJ dml""* \(m n rj 2 ' m^r ' 3 
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FIGURES 

FIG. 1. Derivative of the central component of the NN interaction, rn q ^^- as a function of the 
nucleon-nucleon relative distance in the Skyrme model. Solid curve corresponds to 

m lVrn' With 

U = 93MeV and e = 4.0. Dotted curve represents V^ m " =138) - V^ =0) . 

FIG. 2. Derivative of the tensor component of the NN interaction, rn g §^ as a function of the 
nucleon-nucleon relative distance in the Skyrme model. Solid curve corresponds to w<jf^- with 
f w = 93MeV and e = 4.0. The dashed curve corresponds to the OPEP results. Dotted curve 
represents V ^ =138) - V^ =0) . 

FIG. 3. Derivative of the spin-spin component of the NN interaction, rn q ^^- as a function of 
the nucleon-nucleon relative distance in the Skyrme model. The labels of the curves are the same 
as in Fig. 2. 
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TABLES 

TABLE I. Nucleoli observables (O) and the corresponding derivatives with respect to the quark 



mass. 





U = 54 MeV e = 4.84 


U = 93 MeV e = 4.0 







nl 1 dm q 


O 


m ^ 

1 dm q 


M N 


938 MeV 


52.2 MeV 


1818MeV 


59.6 MeV 


9A 


0.65 


-0.064 


1.01 


-0.069 


gnNN 


11.88 


-0.022 


20.35 


-0.33 



TABLE II. Contributions of the derivative of the NN potential and average kinetic energy to 
the nuclear sigma term. Calculations performed for the Skyrme Model with / 7r =93MeV and e=4.0. 
The upper label in the nucleus symbol denotes the parametrization of the densities according to 
de Jager, de Vries and de Vries[21]. 





E K (MeV) 


Sf (MeV) 


SX(MeV) 


Ty A +K (MeV) 


Ty A +K /A (MeV) 


2 H 


17.8 


-0.86 


-3.60 


-4.46 


-2.23 


3 H 


46.1 


-2.3 


ra-14 


«-16 


ps-5.4 


3 He 


46.1 


-2.3 


«-14 


«-16 


»-5.4 


4 He (3pF) 


106 


-5.09 


-34.5 


-39.6 


-9.9 


10 B (h.o.) 


228 


-10.9 


-107 


-118 


-11.8 


12 C (3 P F) 


274 


-13.1 


-157 


-170 


-14.2 


14 N (h.o.) 


319 


-15.3 


-188 


-203 


-14.5 


16q(3 P F) 


365 


-17.5 


-220 


-238 


-14.9 


20 Ne (2pF) 


456 


-21.8 


-268 


-290 


-14.5 


24 Mg( 3 P F ^ 


547 


-26.2 


-388 


-414 


-17.2 


32g(3pF) 


730 


-35.0 


-557 


-592 


-18.5 


40 Ca (3pG) 


912 


-43.8 


-727 


-771 


-19.3 


56p e (3pG) 


1273 


-61.1 


-1144 


-1205 


-21.5 


118g n (3pG) 


2690 


-129 


-2763 


-2892 


-24.5 


148g m (2pF) 


3374 


-162 


-3556 


-3718 


-25.1 


197 Au (2pF) 


4492 


-216 


-5187 


-5403 


-27.4 


208pk(3pG) 


4742 


-227 


-5274 


-5501 


-26.4 


p(0.17fm- 3 ) 


22.8,4 


-1.09,4 


-38.0,4 


-39.1.4 


-39.1 



14 



TABLE III. Results for the isoscalar parameter of the optical potential including the potential 
sigma term contribution in Eq. (|30|) . The 12 C and 16 experimental values were estimated from 
the analysis described in Ref.[10] and for 40 Ca from the analysis of Ref.[ll]. The value of b^ 1 ^ was 
obtained from the experimental value of the 7r-d scattering length [26] (&0e// = ~~ 0.0099 (m~ 1 )) . 





uV+K 




, (exp) 
°0 r , fr 


2 H 


-0.0025 


-0.0123 


-0.0123(3) 


3 H 


-0.0063 


-0.0162 


- 


3 He 


-0.0063 


-0.0162 


sa-0.015 




-0.0110 


-0.0209 


-0.0212(5) 


10 B (h.o) 


-0 01 SI 


-0 0230 




12 C (3pF) 


-0.0157 


-0.0256 


-0.0260 


14 N (h.o.) 


-0.0160 


-0.0259 




16 (3 P F) 


-0.0165 


-0.0264 


-0.026 (5uMcV, ;-0.027 l29MoVj 


20 Ne (2pF) 


-0.0160 


-0.0259 




24 Mg (3 P F) 


-0.0190 


-0.0289 




32g(3pF) 


-0.0204 


-0.0303 




40 Ca (3 P F) 


-0.0213 


-0.0312 


-0.033 (2bMoVJ 


Pe = 0.17fm" 3 /2 


-0.0216 


-0.0315 


-0.0304 
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